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Director 

Mathematical  and  Information  Sciences 

Air  Force  Office  of  Scientific  Research 

United  States  Air  Force 

Building  410 

Bolling  Air  Force  Base 

Washington,  D.C.  20332 

Attention:  Major  Charles  L.  Nefzger 

AFOSR/NM 

Subject:  Microprocessor  Requirements  for  Implementing  Modern 

Control  Logic,  Contract  F49620-79-C-0078 ,  Technical 
Progress  Report  No.  1  for  the  period  March  1,  1979 
through  August  31,  1979. 

Reference:  Farrar,  F.  A.  and  R.  S.  Eidens:  Microprocessor 

Requirements  for  Implementing  Modern  Control  Logic. 

United  Technologies  Research  Center  Report  R79-944258-2 , 

Final  Technical  Report  prepared  under  Air  Force  Office 
of  Scientific  Research  Contract  F49620-78-C-0017 , 

March  1979. 

Gentlemen: 

la.  The  objective  of  this  research  investigation  is  to  develop 
and  evaluate  analytical  procedures  for  establishing  microprocessor 
requirements  for  implementing  modern  control  logic.  Key  issues  in 
microprocessor  implementation  of  modern  control  logic  include  (1)  accuracy, 
(2)  computat ional  capability,  and  (3)  memory  requirements.  These  require¬ 
ments  must  be  established  for  the  interface  as  well  as  for  the  micro¬ 
processor  control  code. 

lb.  A  Phase  I  first-year  study  directed  toward  establishing 
microprocessor  requirements  for  control  of  linear  systems  was  performed 
at  United  Technologies  Research  Center  (UTRC)  from  1  February  1978  to 

31  January  1979  under  Contract  F49620-78-C-001 7  with  the  Air  Force  Office 
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of  Scientific  Research  (AFOSR) .  This  Phase  1  program  resulted  in  successful 
development  of  analytical  procedures  for  establishing  microprocessor  requirements 
for  multivariable  feedback  control  of  linear  stochastic  dynamic  systems.  The 
developed  procedures  were  evaluated  and  illustrated  by  application  to  (1)  a 
second-order  system  and  (2)  a  linearized  fifth-order  F100  turbofan  engine  model. 
These  results  are  documented  in  the  referenced  report.  based  on  these  results, 
the  subject  contract  was  awarded  to  UTRC  for  a  Phase  II  second-year  program. 

This  second-year  study  is  directed  toward  (1)  validating  Phase  1  results  by 
demonstrating  microprocessor  implementation  of  linear  quadratic  Gaussian  (LQG) 
logic  for  control  of  linear  systems  and  (2)  extending  the  analysis  of  micro¬ 
processor  requirements  for  control  of  nonlinear  systems. 


2.  During  this  report  period  the  UTRC  effort  has  been  directed 

toward  validating  Phase  I  results  by  demonstrating  LQG  control  on  an  Intel 
8080  microprocessor  for  a  second-order  continuous  system.  The  second-order 
system  is  the  same  system  used  in  the  Phase  I  study.  Accuracy,  computational 
capability,  and  memory  requirements  of  the  microprocessor  demonstration 
system  will  be  compared  with  requirements  predicted  using  the  analytical 
procedures  developed  in  the  Phase  1  study.  In  addition,  a  paper  presenting 
the  Phase  1  study  results  was  submitted  and  accepted  for  publication  in  the 
IEEE  Transactions  on  Automatic  Control  and  for  presentation  at  the  18th  IEEE 
Conference  on  Decision  and  Control  (CDC).  Approval  to  submit  the  paper  was 
obtained  from  Major  Charles  L.  Nefzger,  the  AFOSR  contract  monitor.  A  copy 
of  the  cover  letter  and  paper  submitted  for  publication  in  the  Transactions 
is  enclosed  for  Major  Nefzger  with  this  report. 

3a.  To  demonstrate  LQG  control  on  an  Intel  8080  microprocessor  for 

the  second-order  continuous  system  (analog  model),  the  system  shown  in  Fig.  1 
is  being  employed.  The  system  consists  of  the  analog  model,  an  Intel  8080 
microprocessor,  A/D  and  D/A  converters,  and  a  strip  chart  recorder.  The 
continuous  system  dynamics  (analog  model)  are  represented  by 

y+4y  +  2y  =  up;  y(o)  =05 

y(o)  =  0.5  (1) 

zp=  y  +  np 


where  y,  Up,  Zp,  and  Hp  represent  the  system  output,  input,  measurement,  and 
measurement  noise,  respectively.  The  dot  notation  denotes  differentiation  with 
respect  to  time.  The  absolute  maximum  values  are  given  by 


|y|  -  i.o 

|y|  =  1.0 

|y|  =  *o 

|u|max  =  20 


(2) 


The  analog  system  shown  in  Fig.  1  is  scaled  so  that  10  volts  corresponds  to  the 
maximum  value. 
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3b.  For  LQG  control  design,  continuous  system  dynamics  are 

represented  in  state-space  notation.  For  digital  control  implementation, 
system  dynamics  are  normalized  so  that  the  numbers  in  the  control  compu¬ 
tations  range  between  -1.0  and  +1.0.  Normalized  continuous  system 
dynamics  are  described  by  ...... 

X  =  Ax  +  Bu  +  £ 

y  =  Cx  Du  (3) 

2  =  Ex  +  n 

where  the  vectors  x(nxl),  u(mxl),  y(pxl)  and  z(lxl)  represent  the  normalized 
system  states,  inputs,  outputs,  and  measurements,  respectively.  The  random 
process  vectors  e  and  n  represent  white  zero-mean  Gaussian  n-dimens ional 
process  and  1-dimensional  measurement  noise,  respectively.  The  second-order 
system  matrices  (A,  B,  C,  D,  and  E)  are  shown  in  Table  I.  The  companion  form 
is  shown  for  ease  in  relating  the  state  equations  (Eq.  (3))  to  the  system 
input-output  equation  (Eq.  (1)).  (In  the  referenced  report,  the  system 
matrices  were  displayed  in  standard  form.) 

3c.  Normalized  discrete  control  dynamics  implemented  on  the  micro¬ 

processor  are  given  by 


w(k  +  i)  =  </>0d(k)  +  H0z(k  +  i) 
u  (k  +  l)  =  G0  w  (k  +  l )  U) 


The  vector  w  represents  estimated  transformed  state  variables.  The 
matrices  G  and  H  denote  the  deterministic  control  gains  and  filter  gains, 
respectively.  The  gain  matrices  for  the  second-order  system  are  shown  in 
Table  I.  The  transformation  matrix  T  is  an  identity  matrix. 
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4a.  The  software  to  implement  Eq.  (4)  on  a  microprocessor  consists  of 

three  matrix/vector  multiplications  (*^W,  H^Z,  and  G^W)  and  one  vector  addi¬ 
tion  ((♦yW)  +  (HyZ)).  In  addition,  software  for  the  LQG  control  system  includes 
a  routine  to  save  past  state  estimates  and  an  interrupt  service  routine.  The 
interrupt  service  routine  (l)  performs  a  timing  check  to  assure  that  all  con¬ 
trol  computations  are  completed  within  the  sample  time  and  (2)  reads  in  mea¬ 
surement  data  and  updates  the  system  inputs. 

4b.  The  overall  block  diagram  of  the  control  software  described  above 

is  shown  in  Fig.  2.  The  block  diagram  of  the  matrix/vector  multiplication 
code  is  displayed  in  Fig.  3.  Figure  3  indicates  that  several  minor  changes 
(e.g.,  the  order  in  which  pointers  are  initialized  and  updated)  were  made  in 
the  mat r ix/vector  multiplication  block  diagram  presented  in  the  referenced 
report.  These  changes  result  in  more  efficient  microprocessor  implementation. 
Note  that  the  matrix/vector  multiplication  algorithm  is  not  changed  but  rather 
the  way  the  algorithm  is  implemented  has  been  slightly  modified.  Block  dia¬ 
grams  of  the  vector  addition  code,  the  store  state  estimates,  and  the  interrupt 
service  routine  are  shown  in  Fig.  4. 

4c.  Computation  times  as  a  function  of  computer  cycles  to  execute  the 

blocks  of  code  shown  in  Figs.  3  and  4  are  displayed  in  Table  2.  The  cycle  time 
to  execute  a  block  of  code  varies  with  the  microprocessor  used.  Note  also  that 
the  computation  time  is  defined  in  number  of  cycles.  Once  the  number  of  cycles 
for  the  code  on  a  given  microprocessor  has  been  defined  conversion  to  seconds 
for  different  clock  frequencies  is  easily  done  (t(sec)  ■  t (eye les/c lock  frequency 
(hertz) ) . 


5a.  To  compute  the  cycles  for  a  given  microprocessor,  code  must  be 

written  to  execute  the  functions  shown  in  the  block  diagrms  of  Figs.  3  and  4. 
Control  code  for  the  Intel  8080  microprocessor  is  shown  in  Fig.  5.  The  code 
consists  of  software  to  execute  the  control  equations  (Eq.  (4))  and  interface 
software.  A  software  multiplication  subroutine  is  used.  Improvements  were 
made  in  the  preliminary  matrix/vector  multiplication  code  presented  in  the 
referenced  report.  These  improvements  result  in  reduced  execution  time.  The 
code  changes  include  (1)  more  efficient  use  of  the  registers  in  the  multipli¬ 
cation  algorithm  (11. 5X  reduction  in  cycle  time)  and  (2)  more  efficient  memory 
(data  array)  accessing  as  well  as  more  efficient  use  of  the  registers  in  the 
matrix/vector  multiplication  algorithm  (reduction  in  cycle  time  dependent  on 
system  order,  e.g.,  a  24X  reduction  in  the  mat r ix/vector  multiplication  control 
algorithm  is  obtained  for  a  2  x  2  matrix  times  a  2  x  1  vector).  In  addition, 
the  interface  software  was  added  to  the  preliminary  code.  The  code  is  also 
completely  documented.  The  flow  chart  for  the  code  is  shown  in  Fig.  2. 
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Sb.  Table  111  shows  times  associated  with  executing  the  blocks  of  code 

in  Figs.  I  and  .  Using  the  Phase  t  procedure  to  determine  the  computation  time 
for  the  second-order  system  (Tables  11  and  III)  indicates  that  b25J  cycles  are 
required  tor  one  pass  through  the  code.  For  a  2  MU*  clock  (on-board  KOM  memory) 
the  predicted  computation  time  is  equal  to  3.13  ms;  for  a  1.47  MHr  clock  (off- 
board  RAM  memory)  the  predicted  computation  time  is  4.'.'S  ms.  For  convenience 
off-board  RAM  memory  will  be  employed  in  this  study. 

b .  Intel  SObvi  memory  requirements  as  a  function  of  system  state, 

input  and  output  orders  are  shown  in  Table  IV.  These  requirements  include 
the  LQC  control  logic  code  as  well  as  the  interface  logic  code  (interrupt 
service  routine).  Phase  l  procedures  indicat'  that  the  total  memory  required 
for  the  second-order  system  is  •*40  bytes  l  1  S  Pyles  of  RAM  and  -*75  bytes  of 
ROM) . 

7.  Phase  l  study  results  indicate  that  the  Intel  SOSO  with  software 
multiply  can  be  used  to  implement  the  l.q'O  control  law  for  the  second-order 
system.  The  requirements  tor  the  second-order  controller  using  the  standard 
form  are  analytically  computed  to  be  (l)  an  S-bit  word  length  accuracy,  (J)  a 
minimum  sample  time  (based  on  comput  at  tons )  ot  s.J1'  ms  and  a  maximum 

sample  time  (based  on  controller  poles)  ot  700  ms,  and  (3)  -*40  words  of 
memory  (475  words  o t  PROM  and  Is  words  of  RAM).  These  requirements  include 
interface  as  well  as  microprocessor  requirements. 

8.  The  system  shown  m  Fig.  1  has  been  constructed.  The  code  lias 
been  implemented  on  the  Intel  8080  microprocessor.  The  memory  used  agrees 
with  t tie  analytically  predicted  requirements.  The  actual  computation  time 
sampling  interval  has  been  computed  using  an  osc i 1 loscope .  The  minimum  sample 
time  (based  on  actual  code)  was  -i.70  ms.  This  actual  computation  time  is 
approximat e l y  10 »  higher  than  the  analytically  predicted  code.  This  difterence 
is  within  the  ♦  JOT  variation  that  is  expected.  Figure  h  compares  the 

actual  system  response  with  the  Phase- 1  predicted  response.  Figure  b  shows 
that,  as  predicted  in  the  Phase  1  study,  8-bit  accuracy  is  sufficient  for 
implementing  the  second-order  controller. 

4.  In  the  next  report  period  the  second-order  and  fifth-order  controller 

analysis  will  be  completed.  Procedures  for  determining  microprocessor  require¬ 
ments  will  be  extended  to  nonlinear  systems.  In  addition,  results  obtained  in 
the  Phase  l  study  will  be  prepared  tor  present  at  ion  at  the  18th  IEKK  CIV. 
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10.  The  Principal  Investigator,  Mrs.  Florence  Farrar,  is  starting 

full-time  studies  for  her  Ph.D.  at  the  University  of  Connecticut  (UCONN)  in 
September.  Dr.  Robert  Guile,  who  received  his  Ph.D.  from  UCONN  in  1970,  has 
assumed  the  responsibilities  of  Principal  Investigator  effective  September  20, 
1979.  Dr.  Guile  has  most  recently  been  engaged  in  research  in  the  area  of  peak* 
seeking  adaptive  systems.  He  has  previously  been  involved  with  analysis  and 
design  of  digital  electronic  control  systems.  He  has  10  years  of  engineering 
experience,  including  3  years  at  UTRC.  Mr.  James  Krodel  is  continuing  as 
Co-invest igator . 

Very  truly  yours, 

UNITED  TECHNOLOGIES  CORPORATION 
Research  Center 


Assistant  Manager, 
Engineering  Operations 


DPM/dip 
Enc losure 
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TABLE  11 

COMPUTATION  TIME 


Mai r  ix/Vector 

Mult ipl icat  ions 

*> 

n  ♦  ni  ♦  nm 

Funct ion 

Mat  r ix/Vector 

Add  i  t  ion 

n(n-l)  ♦  n(l-l)  ♦  m(n-l) 

Vector  Additions 

n 

Mat  r ix/Vector 

Mult ipl icat ion 

(A  *  )  (N‘  ■*  nt  ♦  mn  -  2n  -m)  ♦  (Ea  >  ♦  La}) 

(2n  *  m)  -3t  ->  ♦  (m  *  L_  , )  (n“  *  nt  ♦  nm) 

♦  ( 2n  ♦  O  ♦  3 

T  ime* 

( Cyc les 1 

Vector  Addition 

State  Store 

(pal  *  Pa2)n_  Pa2  >  pa0 
(Pdl  *  Pd2)n~  Pd2  *  Pd0 

Interrupt  Service 

TC  ♦  (AD) l  ♦  ( DA)m 

*  See  Figs.  3  and  4  for  definitions  of  time  notation  (e.g.,  A  defined  on  Fig.  3). 


R  S'JO-  1 


TABLE  111 


Intel  8080  Execution  Times 
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TABLE  IV 

MEMORY  REQUIREMENTS 


Var i ab  le 

Memory  Type 

Memory  (BYTES) 

A 

Past  state  estimate  (w(k)) 

RAM 

n 

A 

Current  state  estimate  (w(k«-l)) 

RAM 

n 

Measurement  tz(k*l)) 

RAM 

t 

Control  (u(k+l)) 

RAM 

m 

System  matrix  ( $y) 

PROM 

n2 

Kalman  gain  matrix  CHy) 

PROM 

nt 

Control  gain  matrix  (Gy) 

PROM 

mn 

Temporary  Storage 

RAM 

2 

n*"  ♦  3nl  ♦  nm  ♦  mt 

Main  Control  Logic 

PROM 

393 

Computer  Code 

Interrupt  Service 

PROM 

74 

Rout ine 

RAM 

2 

n‘+  3nl+nm  +  mt-*-ra«-  l 

TOTAL 

PROM 

n2  ♦  nl  ♦  mn  ♦  4b 7 

SYSTEM  STRUCTURE  FOR  DEMONSTRATING  INTEL  8080  MICROPROCESSOR  CONTROL 
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INTEL  8080  SOFTWARE  FOR  LOG  CONTROLLER 
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FIG  5  (CONTINUED) 


INTEL  8080  SOFTWARE  FOR  LOG  CONTROLLER 
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FIG  5  (CONTINUED! 


INTEL  6080  SOFTWARE  FOR  LOG  CONTROLLER 
(CONTINUED) 
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INTEL  8080  SOFTWARE  FOR  LOG  CONTROLLER 
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FIG  5  (CONTINUED) 


INTEL  8080  SOFTWARE  FOR  LOG  CONTROLLER 
(CONTINUED) 
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FIG  f.  CONCLUDED* 


INTEL  8080  SOFTWARE  FOR  LOG  CONTROLLER 
(CONCLUDED) 
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